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Abstract 

The potential field approximation has been providing a fast, and computationally 
inexpensive estimation for the solar corona’s global magnetic field geometry for 
several decades. In contrast, more physics-based global magnetohydrodynamic 
(MHD) models have been used for a similar purpose, while being much more 
computationally expensive. Here, we investigate the difference in the field geom¬ 
etry between a global MHD model and the potential field source surface model 
(PFSSM) by tracing individual magnetic field lines in the MHD model from 
the Alfven surface (AS), through the source surface (SS), all the way to the 
field line footpoint, and then back to the source surface in the PFSSM. We also 
compare the flux-tube expansion at two points at the SS and the AS along the 
same radial line. We study the effect of solar cycle variations, the order of the 
potential field harmonic expansion, and different magnetogram sources. We find 
that the flux-tube expansion factor is consistently smaller at the AS than at the 
SS for solar minimum and the fast solar wind, but it is consistently larger for 
solar maximum and the slow solar wind. We use the Wang-Sheeley-Arge (WSA) 
model to calculate the associated wind speed for each field line, and propagate 
these solar-wind speeds to 1 AU. We find a more than five hours deviation in 
the arrival time between the two models for 20 % of the field lines in the solar 
minimum case, and for 40 % of the field lines in the solar maximum case. 

Keywords: Magnetic fields. Models — Solar Wind, Theory — Velocity Fields, 
Solar Wind 


1. Introduction 

The magnetic field above the solar photosphere and in the solar corona is 
weak. Therefore, it is hard to obtain observationally the three-dimensional topol- 


^ Harvard-Smithsonian Center for Astrophysics, 
60 Garden St. Cambridge, MA 02138, 

USA email: ocohen@cfa.harvard.edu 


; p. 1 


SOLA: ms.tex; 3 July 2015 


0:37 





ogy of the magnetic field in the solar atmosphere. Global imaging of the solar 
corona (e.g., TRACE, SOHO EIT, STEREO EUV, HINODE XRT, and SDO 
AIA EUV/X-ray images), clearly shows magnetic loop structure in the corona. 
However, these images indicate the line-of-sight integral of local density and 
temperature structures, so they cannot provide the complete information about 
the magnetic field. 

While limited three-dimensional reconstruction is possible with the extensive 
available solar data, a complete picture of the three-dimensional structure is still 
limited to theoretical models. The most popular method to estimate the global 
three-dimensional coronal magnetic field is the so-called ’’Potential Field Source 
Surface Model” (PFSSM hereafter, Altschuler and Newkirk, 1969 Altschuler 
et ai, 19771. The underlaying assumption in this method is that the large-scale 


structure of the solar corona does not change much over one solar rotation (the 
time-scale for the appearance of new active regions is longer), and it definitely 
does not change within the time it takes perturbations to move along the mag¬ 
netic field (the Alfven time-scale). With that assumption in mind, one can also 
assume that there are no forces in the form of electric currents applyed on the 
magnetic field and that the field is static. Other method to calculate more local 
regions in the corona allows field-aligned currents in the system. This method 
is called the linear/non linear ’’Force-free” method (see review by Wiegelmann 
and Sakurai, 2012[ with references therein). 

The PFSSM assumption that there are no currents in the system immediately 
implies that: 


V X B = 0, (1) 

so the magnetic field vector, B, can be described as a gradient of a scalar 
potential, 


B = -V^-. (2) 

With the additional requirement that there are no magnetic monopoles, V • B = 
0, we can obtain the Laplace equation for the scalar potential: 

= 0. (3) 


Equation can be solved using two boundary conditions. The first boundary 
condition at the solar surface, r = Rq, can be specified using observations of 
the photospheric radial magnetic field (magnetograms). The second boundary 
condition assumes that the coronal magnetic field becomes purely radial at some 
spherical distance, i?ss, with the sphere of r = R^a being known by the name 
’’Source Surface” (SS hereafter). Once the boundary conditions are obtained, 
the solution for the three-dimensional structure of the coronal magnetic field 
can be obtained using harmonic decomposition and the Associated Legendre 
functions, where the harmonic coefficients in such a solution are computed using 
the observed photospheric field. Using well-known recursive relations for the 
polynomial functions and modern computers, the solution can be obtained al¬ 
most instantaneously (see e.g., Luhmann et ai, 2002). This makes the PFSSM 
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extremely popular to use, in contrast to expensive, physics-based models that 
will be discussed below. 

The assumption that the SS is spherical and the assumption about its location 
are the main drawbacks of the method for the following reasons: First, the 
physical interpretation of the mathematically described SS is that it represents 
the surface at which all field lines are supposedly open. In reality, this surface 
is known as the Alfven Surface (AS hereafter), which represents the collection 
of points along the ensemble of coronal field lines at which the solar-wind speed 
exceeds the local Alfven speed [ua — Bj i/47rpj with B being the local magnetic- 
field strength and p being the local mass density. At this point, the wind’s 
dynamic pressure overcomes the magnetic tension in the magnetic loops. As a re¬ 
sult, the field lines are “open“ to the heliosphere and are dragged by the radially 
expanding solar wind. The shape of the AS is far from spherical, and it depends 
on the complexity of the solar magnetic field. Even during solar minimum, where 
the solar magnetic field is nearly dipolar, the AS is not spherical, with the surface 
being much closer to the Sun near the equator (where the magnetic field is weak) 
than above the poles (where the magnetic field is stronger). Despite this trivial 
limitation of the method, very little analytical work has been done to estimate 
the non-hericality of the source surface in the context of the heliosphere (Schulz, 

[199^. 


Frazier, and Boucher, 1978 Levine, Schulz, and Frazier, 1982 Schulz 


The other problem with the outer boundary condition in the PFSSM is the 
assumption where to locate the SS. Traditionally, i?ss is set at 2.5rQ. This 
distance originated from comparison between the height of the coronal loops 
in the PFSSM solutions and in solar eclipse images (Altschuler and Newkirk, 


1969 Hoeksema, 1984). While this value seems to work reasonably well overall, 


it is more than likely that the SS can be located further out, in particular when 
the magnitude of the photospheric field increases. Recently, DeForest, Howard, 
and McComas (2014) have used STEREO-A/COR2 data to identify returning 
structures in the solar corona and solar wind. These returning structures are 
indicators for sub-Alfvenic regions, meaning that these regions are below the AS. 
They found that overall, the AS is located at least above Gr©, and it may even 
be located above 12rQ for coronal holes and above ISr© for helmet-streamers. 

The PFSSM is commonly used as the initial condition for magnetohydro¬ 
dynamic (MHD), physics-based models. The main advantage of these models 
over the PFSSM model is that they provide solutions for all the MHD variables 
of density, velocity field, magnetic field, and pressure. However, as mentioned 
above, they are much more computationally expensive as they require a full 
three-dimensional computation for all the variables. Riley et al. (2006) performed 
a detailed comparison between MHD and PFSSM solutions for the global solar 
corona. They used a simplified, MHD model with a simple adiabatic energy equa¬ 
tion (constant polytropic index of F = 1.05, and no sources: Mikic and Linker, 
1994[ Lionello, Mikic, and Schnack, 1998 Linker et al, 1999 Mikic et al., 1999) 


to compare solutions driven by magetograms taken during both solar minimum 
and solar maximum periods. They found that the coronal holes structure in the 
MHD and PFSSM solutions were tentatively similar, while the coronal holes’ 
area was consistently smaller in the PFSSM solutions. This suggests that some 
of the magnetic flux that is being opened by the solar wind in the MHD model 
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remains closed in the PFSSM model. One can argue that the simplified MHD 
model for the solar wind may not realistically energize the solar wind and that 
with a more realistic model developed by the same group, such as Lionello, 
Linker, and Mikic (2009) or Riley et al. (2011), the discrepancies in coronal hole 
mapping may be even greater. Unfortunately, the MHD solution for the solar 
wind itself was not presented by Riley et al. (20061. 

Riley et al. (2006) have also compared the shape of the SS with the shape 
of an iso-surface extracted from the MHD solution for nearly radial field (i.e. 
|i?r|/|S| = 0.97). They found that the solar minimum configuration is consis¬ 
tent with analytical predictions of non-spherical PFSSM (Schulz, Frazier, and 
Boucher, 1978[ Levine, Schulz, and Frazier, [l982 Schulz, [l997 ), while the solar 
maximum MHD iso-surface is essentially spherical when averaging the effect of 
local, time-dependent active regions. We would like to stress that this radial- 
field criteria may be misleading in the physical context of the SS. As mentioned 
above, we need to consider the point where the field line is open, not where the 
field is radial, as a radial field can still be closed. For example, the top of the first 
open helmet-streamer field line can become radial at a relatively close distance 
from the Sun, but it is open by the solar wind at a greater distance. Indeed, the 
near radial sphere found by Riley et al. ( 2006) is at much greater distance than 
the Rss = 2.5rQ in the PFSSM. 

Lee et al. (2011) compared the observed magnetic flux at 1 AU with the 
magnetic flux predicted by the PFSSM using different locations of the SS. They 
found that in order to better match the observed and predicted fluxes, the SS 
should be placed below 2.5rQ (1.5 — 1.8rQ). This is clearly due to the fact that 
the PFSSM is a static model that truncates the loops at certain height so that 
the lower the SS is, the more open flux is added. Some MHD models (e.g. Cohen 
et al., 2007), truncating the magnetic loops below the height at which they 
are actually opened by the solar wind, may over-energize the wind, leading to 
unrealistic solar-wind speeds in the solution. This is due to the fact that the solar- 
wind energization is specified according to the initial, potential-field distribution. 

An important application of the PFSSM has been introduced by Wang and 
Sheeley (1990), where they presented empirical inverse relation between the 
terminal solar-wind speed, Usw, and amount of expansion of the magnetic-flux 
tube from which that particular wind originated. The expansion factor, /g, is 
defined as the ratio between the magnetic flux at the flux tube’s footpoint to 
that at the SS: 


_ B{rQ)R% 

Js — 


B{Rss)Rl 


( 4 ) 


SS 


for a given magnetic-field line. The value of fs is calculated from the PFSSM 
solution. By fitting their predicted solar-wind speed to solar wind data at 1 AU, 
Wang and Sheeley (1990) derived the following expression for the solar-wind 


speed as a function of fs'. 


= u„ 




-1 


( 5 ) 


Here Umin is the minimum observed solar-wind speed, A is a coefficient, and I 
is a power-law index. Since its initial release. Equation has been improved to 
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include real-time updating of the input magnetogram (Arge and Pizzo, 20001 
and a smoothing factor that depends on the angular distance of the particular 
field line from the nearest coronal hole boundary, Oi,. With this new factor and 
the real-time update, Equation takes the following form: 


= Un 


■ A(1 ■ 




exp[l- (6»b/D)“]}", 


( 6 ) 


where the different parameters have been continuously changed over the years 
to improve the fitting to solar-wind data. Equation is known as the Wang- 
Sheeley-Arge (WSA) model which predicts the solar-wind speed at 1 AU at a 
given time from magnetogram data and the PFSSM (McGregor et ai, 20111. 

The WSA model is widely used by the heliophysics community. It has been 
used for years at NOAA’s Space Weather Prediction Center (SWPC) and it 
has also been used to specify the solar-wind speed at the inner boundary of 
the ENLIL code (Odstrcil, [2003 ): an MHD code that propagates the solar wind 
from the solar corona to 1 AU. The WSA-ENLIL combination has been used 
for many studies of the solar wind directly by individual researchers and via the 
Community Coordinated Modeling Center (CCMC). Despite its extensive usage, 
the WSA model depends heavily on the PFSSM and the determination of /g. It 
is also known to be sensitive to the magnetogram observatory and resolution as 


a result of different mapping of the potential field (e.g. Lee et al, 2009 2011). 
Recently, Poduval and Zhao (2014) have studied the different in the expected 


fs between PFSSM and the current sheet source surface (CSSS) model, which 
analytically allows field lines to change their expansion shape between tq and 
Rss ■ They found that this geometrical correction improves the prediction of the 
solar-wind speed at 1 AU and that the WSA-ENLIL approach can be improved 
beyond the simple PFSSM. 

In this article, we study the deviations between the PFSSM prediction for fs, 
and fs calculated at the Alfven point from the MHD solution for individual field 
lines. The MHD model accounts for the acceleration of the solar wind and the 
heating of the corona self-consistently. Therefore, it captures major deviations in 
the PFSSM field expansion and demonstrates whether or not the PFSSM wind 
prediction is quantitatively reasonable. 

In the next section, we describe the MHD model and data that we use to 
obtain fs- We show the results in Section and discuss their implications in 
SectionWe conclude our findings in Section]^ 


2. Model and Procedure 


The solar-wind MHD solution is obtained using the BATS-R-US generic MHD 

Toth et al., 


code (Powell et al., 1999 Toth et al, 2005 


2012). The particular 


version of the code for simulating the solar corona and the solar wind includes 
physics-based coronal heating and acceleration of the solar wind by dissipation of 
the turbulent energy of Alfven waves, as well as thermodynamics effects such as 
radiative cooling and electron heat conduction (Oran et al., 2013 Sokolov et al.. 


2013 van der Holst et ai, 2014). The initial condition for the three-dimensional 
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magnetic field is obtained using PFSSM extrapolation of magnetogram data as 
described above (with the SS set at 2.5rQ). The initial potential field becomes 
non-potential as the corona is heated by the dissipation of turbulent energy of 
Alfven waves, and the solar wind is accelerated by the increased thermal and 
wave pressure gradient. The steady-state solution includes a hot corona and a 
fully developed solar wind (e.g. Pneuman and Kopp, 1971 Rusin et al., 2010), 
which stretches the initial magnetic field into helmet-streamer topology, which 
resembles the coronal structure observed during solar eclipses. 

The model uses a spherical grid that is stretched in the radial direction. 
The smallest grid size (at the inner boundary) is O.OItq in the radial direction 
and O.O 2 r 0 in the angular direction, and it increases to 0.5rQ near the outer 
boundary of 24rQ. We implement adaptive mesh refinement (AMR) near the 
coronal current sheet as its location dynamically changes by the wind stretching 
the field lines; the resolution around the current sheet is O.lrQ. A convergence is 
typically achieved after 60,000 iterations, where this rather long time is mostly 
due to the detailed, second-order radiation terms that are included in the model. 
All of the simulations were performed on NASA’s PLEIADES supercomputer at 
NASA AMES using 256 processors. With this setting, each solution is obtained 
within about six hours of physical time. 

Once a steady-state is obtained for a particular case (see Section]^ below), we 
identify all of the points in the solution that belong to the AS, with a selection 
condition that the Alfven Mach number Ma = Usvj/ua = 1 ± 0.001. Typically, 
we identify between 1000 to 2000 points, depending on the shape of the AS 
and the grid structure on it. Eor each identified point on the AS, we begin an 
Euler-based field-line tracing procedure towards the solar surface. The field-line 
tracing defines a spatial path step [±dsi] in the direction of the unit magnetic 
field vector [bp] for the particular point [ri]. The magnitude of dsi is the average 
distance between the surrounding cell centers, and the sign is determined so that 
the tracing is done towards, not away from the solar surface. We then map the 
end location of r' = ri -|- dsibi into the center of the nearest grid cell, ri+i. The 
next step is then calculated using r' = ri+i -|- dsi+ibi+i and so forth until the 
tracing crosses the solar surface. During the field-line tracing procedure, we save 
the following information: i) the location of the starting point on the AS, the 
radial distance of the starting point from the Sun, and the field magnitude at 
that point; and ii) the location and field magnitude at the field line footpoint on 
the solar surface. Erom each identified field-line footpoint, we begin a field-line 
procedure in the PFSSM solution until we reach the SS located at r = 2.5rQ 
and we save the field magnitude at that point. All of these parameters enable 
us to calculate the expansion factor at the SS [fss], and the expansion factor at 
the AS [/as] (with the Alfven point radius [i?as] and the field magnitude at that 
point [Ras] replacing the radius and the field strength of the SS in Equation]^. 

Once the field-line tracing procedure is complete, we have all of the infor¬ 
mation needed to compare the flux-tube expansion between the SS and the 
AS. 
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3. Selected Data 


As the WSA model is known to be sensitive to different magnetogram input, we 
investigate three effects on the expansion factor differences between the SS and 
the AS: First, we investigate the differences between solar-minimum and solar- 
maximum periods using MDI data (www.sun.stanford.edu) for Carrington Rota¬ 


tion (CR) 1922 (for ranges April-May 1997, solar minimum) and CR 1958 (for 
ranges January 2000, solar maximum). Second, we study the effect of the order of 
the harmonic expansion [n] on the differences in the expansion factor, where we 
use the two MDI magnetograms decomposed with n = 90 and n = 180. Third, we 
compare four different magnetograms datasets for CR 2104 (October-December 
2010, rising phase) observed by SDO HMI with high (3600x1400) and low 
(720x360) resolution (http://hmi.stanford.edu/data/synoptic.html |, MDI, and SO-| 
LIS (http://solis.nso.edu). These four datasets are decomposed with n = 180. 

It should be mentioned that, overall, there is a flux difference of about 20 % 
between the MDI and HMI low-resolution magnetogram, and 10 % between MDI 
and the HMI high-resolution map. The flux of the SOLIS map is significantly 
lower than all the other maps. Overall, all of the magnetic features in the SOLIS 
map are almost 80 % weaker than in the MDI map. 


4. Results 


Figure demonstrates the similarity of the three-dimensional field geometry 
between the PFSSM and the MHD model below the SS. Each pair of three- 
dimensional field lines are selected at the same footpoint, and are then drawn 
in each of the three-dimensional solutions. It can be seen that the geometry of 
most of the field lines, which are not associated with closed loops is very similar 
in both models below the SS. The main notable differences are for the field lines 
that belong to closed loops or overlaying closed loops. This is due to the fact 
that the solar wind existing in the MHD solution stretches and opens some of 
the closed loops in the PFSSM. 

Figure shows meridional slices of the MHD solutions using MDI magne¬ 
tograms for CR 1922 and CR 1958, and an HMI high-resolution magnetogram 
for CR 2104. The slices are colored with solar-wind speed contours, along with 
the AS indicated as black line. The dipolar topology of solar minimum period 
(CR 1922) is seen in the solar-wind structure, along with a thin, equatorial region 
of slow solar wind near the current sheet, and a rather symmetric AS. During the 
solar maximum period (CR 1958), the field topology is more complex, leading 
to an extensive slow-wind region and a tilted current sheet. The field complexity 
is also reflected in the shape of the AS. CR 2104 (November-december 2010) 
occurred at the beginning of the rising phase after the extended solar minima at 
the end of Solar Cycle 23. The dipolar topology and the aligned current sheet 
are still present in the solution, while the region occupied by slow wind is more 
extended than in the solution for CR 1922. The shape of the AS represents an 
intermediate state between solar minimum and solar maximum. For all solutions, 
the distance of the AS is in an overall agreement with DeForest, Howard, and 


McComas (2014). 
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Figure 1. Selected three-dimensional field lines, extracted from the same footpoint for the 
PFSSM (blue) and MHD solution (red) for CR 1958 MDI magnetogram with n=180. The inner 
sphere represents the solar surface and the solid-white line represents the SS at r = 2.5rQ. 


Figure shows how the flux-tube expansion changes from the solar surface to 
the AS for two selected flux tubes. One flux tube (left) gradually expands with 
rather constant rate. This topology is typical of coronal holes and fast-wind 
regions. The flux tube on the right represents open flux tube at the boundary 
of a closed loop, which is a typical source region for the slow solar wind. Here 
the expansion significantly changes from the lower regions to the AS. Figure]^ 
demonstrates the physical essence of the relation between the solar-wind speed 
and fs shown in Equation and Equation 

4.1. Tracing the Field Lines in the MHD Model and the PFSSM 
from the Same Footpoint 

The results of the field-line tracing are shown in Figure and Figure Each 
point in the top four scatter plots represents the value of the expansion factor 
at the SS [fss] vs. the value of the expansion factor at the AS [fas]- Figure]^ 
shows the results for MDI magnetogram data taken during CR 1922 and CR 
1958 with n=90 and n=180. Figure shows the results for magnetogram data 
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Table 1. Linear fit and correlation parameters for field line tracing from the SS to the AS 


Magnetogram 

Input 

a 

P 

h 

Correlation 

7 

S 

u 

Correlation 

MDI CR 1922 n=90 

12±2 

0.7±0.19 

0.1 

25±24 

0.89±0.04 

0.6 

MDI CR 1922 n=180 

5±2 

1.6±0.20 

0.3 

220±17 

0.89±0.03 

0.7 

MDI CR 1958 n^90 

95±5 

0.2±0.01 

0.4 

143±10 

0.71±0.02 

0.6 

MDI CR 1958 n^l80 

116±6 

0.1±0.20 

0.20 

184±13 

0.63±0.03 

0.4 

HMI high-res CR 2104 

51±3 

0.1±0.01 

0.3 

199±10 

0.58±0.02 

0.6 

HMI low-res CR 2104 

35±4 

0.2±0.02 

0.4 

122±13 

0.78±0.03 

0.8 

MDI CR 2104 

35±5 

0.4±0.02 

0.5 

38±9 

0.93±0.02 

0.9 

SOLIS CR 2104 

37±2 

0.3±0.02 

0.4 

136±11 

0.71±0.02 

0.7 


taken by the four magnetogram sources for CR 2104. The bottom four scatter 
plots in each figure show the wind speed associated with a given value of the 
expansion factor. The associated speed is calculated using the relation: 


240 + ' 


500 


(l + /s)0.22_ kms \ (7) 

There are a number of different coefficients in the literature that were modified 
to tune the model prediction of the solar-wind speed at lAU. Here, we choose 
the av erage coefficients betwee n the different published values (e.g., Arge and 
Pizzo, 2000 McGregor et ai, 2011 the value of A = 500 km s~^ is also a 


representative mid-point between the fast and the slow solar wind). In order to 
isolate the effect of the expansion factor, we also neglect the smoothing factor 
in Equation that depends on the distance from the coronal hole boundary as 
it is important mostly for real-time tracking of the solar wind from the Sun to 
the Earth. 

For each scatter plot in Figure and Figure we calculate a linear fit in the 
form of fas = Oi + /3/ss, and Uas = 7 + Suss- We also calculate the errors of these 
coefficients, as well as the correlation parameters, and all of these parameters are 
summarized in TableOverall, the correlation for the expansion factors ranges 
between 0.1-0.4, with the correlation being notably smaller for the MDI case of 
CR 1922 with n=90. The correlation for u ranges between 0.44-0.86. The errors 
of the slope /3 for the MDI cases are higher than the cases for CR 2104, with the 
exception of CR 1958 with n=90. 

Overall, the the correlation for the expansion factor between the two models 
is poor and the scatter of the plots is quite large, even for the ’’simple” solar 
minimum configuration. The correlation for the associated wind speeds is better, 
but it is still less than 0.7 for most of the magnetogram sources. Since the 
expansion factor represents the field geometry, this low correlation is another 
indication that the three-dimensional field geometry in the MHD solution is 
significantly different than that of the PFSSM. The parameters presented in 
Table [^provide scaling parameters between the two models for selected cases of 
magnetogram sources. More extensive calculations of longer periods and for all 
magnetogram sources are left for future work. 
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Table 2. Percentage of the magnitudes of /as vs. 
fss for field line tracing from the SS to the AS 


Magnetogram Input 

/as > fss 

fss ^ /as 

MDI CR 1922 n=90 

95% 

5% 

MDI CR 1922 n=180 

97.5% 

2.5% 

MDI CR 1958 n=90 

36% 

64% 

MDI CR 1958 n=180 

40% 

60% 

HMI high-res CR 2104 

35% 

65% 

HMI low-res CR 2104 

35% 

65% 

MDI CR 2104 

18% 

82% 

SOLIS CR 2104 

50% 

50% 


Table shows the percentage of field lines for which /as is greater than f^s 
and vica versa for all cases. Overall, /as is greater than fss for solar-minimum 
periods, and fss is greater than /as for solar-maximum periods. The latter trend 
appears also for CR 2104 (rising phase of Solar Cycle 24). 

4.2. Assuming a Purely Radial Field Above the Source Surface 

The WSA model commonly produces a spherical distribution of the expan¬ 
sion factor and the associated wind speed at the SS (i.e., fs{Rss,(^,4’) 

4’))- The distribution of Usw is then propagated ballistically to predict 
the wind speed distribution at 1 AU, Usw{r = lAU, 9, (f)), assuming that the field 
is purely radial above the SS. We can take advantage of the data calculated 
here to test this approach and compare it with the tracing of field lines in both 
models from the same footpoint. Here, instead of starting our tracing at the AS, 
we start at a given point on the SS, and trace it back to the Sun in order to find 
its footpoint. This tracing identifies ’’open” field lines in the PFSSM. Instead 
of tracing the field lines between the MHD model and the PFSSM model, we 
now simply find the point on the AS with the same angle as the point on the SS 
((^’asj '('i) and {rss,9i,4>i)). Once the three points are identified, we calculate 

the expansion factor in the same manner as in the previous section. 

Figure and Figure show the result for the this radial propagation from 
the SS to the AS. A notable feature is that these results are less scattered than 
the results for field line tracing between the models. Tables and show the 
properties of the linear fit and correlation functions for the radial propagation. 
Overall, the correlation between the values of fs at the SS and the AS are better 
than the full field tracing, and almost all the correlations for fs are about 0.5 or 
higher (with the exception of CR 1958 MDI with n=180), better than the values 
presented in Tablej^ Tablej^shows the percentage of cases at which /as is larger 
than fss and visa versa. It is consistent with the values appear in Table except 
that here, fss > /as in almost all the points in the cases for CR 2104. 

The result of these calculations could be used to get better predictions from 
the common WSA approach based on PFSSM only. 
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Table 3. Linear fit and correlation parameters for radial propagation from the SS to the AS 


Magnetogram 

Input 

a 

p 

h 

Correlation 

7 

6 

u 

Correlation 

MDI CR 1922 n-90 

2.1±0.2 

0.94±0.03 

0.76 

-lOitll 

1.00±0.02 

0.87 

MDI CR 1922 n-180 

2.8±0.3 

0.83±0.03 

0.64 

181±17 

0.65±0.03 

0.59 

MDI CR 1958 n-90 

48.9±2.6 

0.13±0.01 

0.48 

141±8 

0.75±0.02 

0.69 

MDI CR 1958 n-180 

63.4±2.3 

O.lOitO.Ol 

0.29 

202±12 

0.62±0.03 

0.46 

HMI high-res CR 2104 

1.0±0.1 

0.82±0.01 

0.53 

69.8±4.7 

0.88±0.01 

0.80 

HMI low-res CR 2104 

20.Oil.7 

0.18±0.01 

0.65 

56.7±7.8 

0.95±0.02 

0.94 

MDI CR 2104 

2.0±0.1 

0.55±0.01 

0.91 

212.5±8.9 

0.63±0.02 

0.94 

SOLIS CR 2104 

19.0±1.0 

0.19±0.01 

0.63 

101.2±6.6 

0.83±0.12 

0.89 


Table 4. Percentage of the magnitudes of /as vs. 
fss for radial propagation from the SS to the AS 


Magnetogram Input 

/as ^ fss 

fss ^ /as 

MDI CR 1922 n=90 

82% 

18% 

MDI CR 1922 n=180 

77% 

23% 

MDI CR 1958 n=90 

13% 

87% 

MDI CR 1958 n=180 

25% 

75% 

HMI high-res CR 2104 

6% 

94% 

HMI low-res CR 2104 

1% 

99% 

MDI CR 2104 

1% 

99% 

SOLIS CR 2104 

5% 

95% 


5. Discussion 

The ultimate goal of the investigation presented here is to quantify the change 
in the predicted solar-wind speed as a function of the expansion factor when 
the latter is calculated at the SS and at the AS. Our results show consistent 
deviation between fss and /as, ranging from 10-50% difference and even more. 


5.1. The Effect of Solar Cycle Variations and the Order of the 
Harmonic Expansion 

Figure shows that the the results for CR 1922 (solar minimum) represent 
mostly field lines with very small expansion and associated fast solar-wind speed 
above 500 km s~^, and a rather small scatter around the linear-fit line. This is 
a good representation of the rather stable observed fast solar wind (McComas 


et al., 2007). For the vast majority of these field lines (i.e. for the fast solar wind), 


the PFSSM under-estimates the expansion factor at the Alfven surface. This is 
due to the fact that while the expansion is small in these nearly parallel flux 
tubes, their geometry is changed by the solar wind beyond the source surface 
and some further expansion occurs. 

In contrast, the results for CR 1958 (solar maximum) show that the major¬ 
ity of the held lines have large expansion factor (above 100) , slow associated 
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solar-wind speed, and a large scatter around the linear-fit line. This is a good 
representation of the sporadic observed slow solar wind (McComas et ai, 2007). 
The majority of these field lines show fss greater than /as, which means that 
for the slow solar wind, the PFSSM over-estimates the expansion factor. This is 
the result of the artificial truncation and opening of the closed loops at the SS, 
which forces the field lines to be radial at the SS. In the more realistic, MHD 
solution, these field lines are not necessarily radial at the SS, and the result is 
an expansion factor that is consistently larger in the PFSSM comparing to the 
MHD model for large values of /g. 

Interestingly, CR 2104, which represents the rising phase, shows a trend sim¬ 
ilar to the solar-maximum case, as it already introduces significant amount of 
field lines associated with slow solar wind. Both of the trends mentioned above 
provide an indication that the solar wind continues to accelerate and reshape 
the coronal magnetic field beyond the SS. It has been suggested by Gilbert, 
Zurbuchen, and Fisk (2007) that instead of imposing the SS at a particular 
location, the field needs to relax to magnetic-pressure equilibrium in order to be 
fully radial. 

All the trends mentioned above also appear in the comparison between the 
values at the SS and AS for the same angle (radial propagation). Here as well, 
fs{SS) over-estimates fs{AS) for solar maximum and the rising phase. From 
these results we conclude that the expansion factor is similar at the SS and the 
AS, only during pure solar minimum conditions, when a significant amount of 
solar wind comes from coronal holes (fast wind). The flux tubes associated with 
this kind of wind are nearly parallel to each other, and the assumption of a pure 
radial field that drops with (/g « 1) is valid. For other times, when most of 
the wind comes from the regions of streamers/pseudo-streamers, the field is far 
from radial above the SS. 

It seems like the radial propagation smears up the field lines at which the 
non-radial drop in the field is the greatest. As a result, the scatter around the 
linear fit is smaller, and the correlation is better than the full field line tracing. 

The order of the spherical-harmonic expansion shows very little difference in 
the overall MDI solutions. Both the linear fit and the comparison between the 
magnitudes of fss and /as are quantitatively similar. This means that despite 
the difference in the field geometry due to the different resolution, the overall 
mapping of the field lines in the PFSSM comparied to the MHD model is consis¬ 
tent for each line, regardless of the number of field lines (i.,e. the magnetogram 
resolution). 


5.2. The Effect of Magnetogram Source 

The solutions for CR 2104 show very little difference between the magnetogram 
sources. This is quite consistent with the result above that the solution is not 
sensitive to the order of the harmonic expansion. The solutions for the HMI 
magneto grams are almost identical, while the MDI and the SOLIS solution 
show similar trends, but with slightly different parameters. The most notable 
difference appears in the comparison between the magnitude of fss and /as for 
the SOLIS magnetogram solution, where one is equally bigger than the other. 
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Table 5. Fraction of field lines with solar-wind speed ratio above 
10 % and 20 % for field line tracing from the SS to the AS 


Magnetogram 

Input 

solar-wind speed 
Ratio Above 10 % 

solar-wind speed 
Above 20 % 

MDI CR 1922 n=90 

18% 

6% 

MDI CR 1922 n=180 

20% 

5% 

MDI CR 1958 n=90 

35% 

16% 

MDI CR 1958 n=180 

43% 

26% 

HMI high-res CR 2104 

26% 

9% 

HMI low-res CR 2104 

26% 

6% 

MDI CR 2104 

19% 

4% 

SOLIS CR 2104 

20% 

6% 


This is probably due to the fact that the overall magnetic flux and the magnitude 
of the strong-held features is lower than the other magnetograms. 

5.3. Impact on the Solar Wind Arrival Time to Earth 

Beyond its physical interpretation, the relation between the expansion factor 
and the solar-wind speed is used by the WSA model to predict the arrival time 
of a solar wind element to the Earth. We use simple ballistic propagation of the 
solar-wind speed obtained from Equation to calculate the solar-wind travel 
time from the SS, and from the AS to 1 AU. Figure shows the arrival-time 
difference (in hours) as a function of the ratio Uss/uas for a given magnetic held 
line. Most of the points are located between 10 % or less speed difference (ratio 
between 0.9 to 1.1), and an arrival-time difference of about hve hours or less. 
However, a signihcant fraction of points have speed difference of 20 % or more 
(ratio of less than 0.8 or more than 1.2), and an arrival-time difference of more 
than ten hours. The fractions of points with a speed ratio above 10 % and above 
20 % for all cases are summarized in Teble It shows that for solar minimum 
and fast-solar-wind speeds, the arrival time difference is more than hve hours 
for about 20 % of the held lines, and for about 40 % (nearly half) of the held 
lines for solar maximum and slow solar-wind speeds. The solutions for CR 2104 
seems to represent an intermittent stage. 

Figurej^shows a similar calculation for the radial propagation from the SS to 
the AS. Here, the difference in the arrival time for CR 1922 becomes very small 
comparing to these cases with the full held line tracing. For the other cases, 
the solar-wind speed ratio appears mostly for the values of less than one, in 
contrast to the symmetric structure in Figure This is due to the consistency 

of MSS) > MAS). 

6. Summary and Conclusions 

In this article, we investigate the difference in the held geometry between a global 
MHD model and the potential held source surface model (PFSSM) by tracing 
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Table 6. Fraction of field lines with solar-wind speed ratio above 
10 % and 20 % for radial propagation from the SS to the AS 


Magnetogram 

Input 

solar-wind speed 
Ratio Above 10 % 

solar-wind speed 
Above 20 % 

MDI CR 1922 n=90 

2% 

1% 

MDI CR 1922 n=180 

45% 

2% 

MDI CR 1958 n=90 

9% 

2% 

MDI CR 1958 n=180 

46% 

27% 

HMI high-res CR 2104 

39% 

15% 

HMI low-res CR 2104 

48% 

12% 

MDI CR 2104 

41% 

6% 

SOLIS CR 2104 

41% 

3% 


individual magnetic field lines from the Alfven surface (AS), through the source 
surface (SS), all the way to the field line footpoint on the solar surface, and then 
trace them back to the SS in the PFSSM. We also compare two points along 
the same radial line at the SS and the AS. We study the effect of solar-cycle 
variations, the order of the harmonic expansion, and the magnetogram source. 

We find that for solar-minimum periods, the solutions are associated with 
small values of fs and fast solar-wind speeds, and the PFSSM under estimates 
the flux-tube expansion due to the further change of the filed geometry beyond 
the SS. For the solar-maximum case, the solutions are associated with large 
values of fs and slow solar-wind speeds. Here, we find that the PFSSM over¬ 
estimates the expansion factor due to the truncation of the magnetic loops and 
the radial forcing of the field at the SS, which leads to a consistent change in the 
field geometry between the potential field solution and the steady-state MHD 
solution. 

Finally, we estimate the difference in the solar-wind arrival time to 1 AU 
between the two models. We find that the deviation of the arrival-time between 
the two models for 20% of the field lines in the solar-minimum case, and for 
40 % of the field lines in the solar-maximum case is more than five hours. 
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Figure 2. Three meridional slices colored with solar-wind speed contours for MDI CR 1922 
(left), MDI CR 1958 (middle), and HMI high-resolution CR 2104 (right). Selected three-dimen¬ 
sional field lines (not two-dimensional projection) are shown in white, and the AS is shown as 
solid-black line. The three yellow circles represent distances of 2.5, 10, and 15r0, respectively. 
The size of the displayed box is 17.5r0 in each direction. 
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Figure 3. Two selected flux tubes with their boundaries shown in yellow, and their expansion 
from the solar surface to the AS. The inner sphere represents the solar surface, the outer, 
transparent sphere represents the SS at r = 2.5r0 (both spheres are colored with the radial 
field strength), and the blue ellipse represents the area of the flux tube at the AS. 
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Figure 4. Scatter plots showing the value of the expansion factor at the SS and at the AS, 
based on field line tracing, for MDI magnetograms taken during CR 1922 and CR 1958 using 
n=90 and n=180. Also shown is a comparison between the wind speed associated with each 
value of the expansion factor based on Equation]^ Each panel includes a linear fit marked as 
solid red line. The parameters of the linear fit for each case are summarized in Table ^ 
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Figure 5. Scatter plots showing the value of the expansion factor at the SS and at the 
AS, based on field line tracing, for magnetograms taken during CR 2104. Also shown is a 
comparison between the wind speed associated with each value of the expansion factor based 
on Equation[^ Each panel includes a linear fit marked as solid red line. The parameters of the 
linear fit for each case are summarized in Tabled 
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Figure 6. Scatter plots showing the value of the expansion factor at the SS and at the AS, 
based on radial propagation, for MDI magnetograms taken during CR 1922 and CR 1958 using 
n=90 and n=180. Also shown is a comparison between the wind speed associated with each 
value of the expansion factor based on Equation]^ Each panel includes a linear fit marked as 
solid red line. The parameters of the linear fit for each case are summarized in Table 
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Figure 7. Scatter plots showing the value of the expansion factor at the SS and at the 
AS, based on radial propagation, for magnetograms taken during CR 2104. Also shown is a 
comparison between the wind speed associated with each value of the expansion factor based 
on Equation]^ Each panel includes a linear fit marked as solid red line. The parameters of the 
linear fit for each case are summarized in Table [s] 
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Figure 8. Differences in wind arrival time to Earth as a function of the ratio between the 
wind speed calculated at the SS and at the AS (uss/u.as), based on field line tracing. 
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Figure 9. Differences in wind arrival time to Earth as a function of the ratio between the 
wind speed calculated at the SS and at the AS (ttss/uas), based on radial propagation. 
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